Alignment-to-orientation conversion and nuclear quadrupole resonance 
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The role of alignment-to-orientation conversion (AOC) in nuclear quadrupole resonance (NQR) 
is discussed. AOC is shown to be the mechanism responsible for the appearance of macroscopic 
orientation in a sample originally lacking any global polarization. Parallels are drawn between NQR 
and AOC in atomic physics. 
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The phenomenon of alignment-to- orientation conver- 
sion (AOC) Him SUB] has been of recent interest 
to the atomic physics community because it is an im- 
portant physical mechanism in experiments involving the 
evolution of atomic ground-state polarization^ in exter- 
nal fields. In a simple example of atomic AOC, opti- 
cal pumping by linearly polarized light produces align- 
ment (i.e., the rank-two quadrupole moment, which has 
a preferred axis but no preferred direction) in an ini- 
tially isotropic atomic ground state. Application of a 
static electric field along a direction other than that of 
the atomic alignment axis induces quantum beats that 
result in orientation (i.e., the rank-one dipole moment, 
which has angular momentum biased in one direction). 

This Letter draws an analogy between the quadratic 
Stark splitting encountered in atomic physics experi- 
ments and the nuclear quadrupolar coupling encountered 
in solid-state nuclear quadrupolar resonance (NQR) 0,13, 
and nuclear magnetic resonance (NMR) experiments. 
In particular, we show that in a pulsed NQR experiment, 
bulk magnetization is created via AOC in an ensemble 
initially having zero net polarization. We present an an- 
alytic calculation of the NQR signal produced by a pow- 
der consisting of randomly oriented crystallites and use 
techniques to visualize nuclear polarization that were de- 
veloped to aid in the understanding of polarized atomic 
systems. Given the similarities between atomic physics 
and NQR/NMR experiments, other opportunities may 
exist to use knowledge of one of these two well-developed 
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^ Here, we use polarization as a generic term describing the 
anisotropy of a quantum system such as a nucleus or atom. Var- 
ious types of polarization can be described by polarization mo- 
ments associated with corresponding spherical tensor ranks. 

^ Studies of NQR have intensified in recent years, with NQR in 
such 7 = 1 nuclei as ^ "^N and finding applications in, for 
example, biochemis try LlOl , and i n e xplosives, land mine, and 
narcotics detection [TuliaiiallTHTsl . 



fields to provide insight into the other, or to aid in the 
design of new experiments. For example, recent work in 
atomic physics on the selective creation and detection of 
various multipole moments |l(i| could be adapted for use 
in nuclear systems. 

We first describe the relationship between the Stark 
and nuclear quadrupole Hamiltonians, and then go on 
to discuss NQR dynamics in more detail. The effect of 
an electric field E, directed along the quantization axis, 
on atoms of angular momentum F is described by the 
familiar quadratic Stark Hamiltonian: 
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Here aa and a2 are the scalar and tensor electric polar- 
izabilities of the atoms, respectively. The scalar polar- 
izability term causes a uniform shift of all the magnetic 
sublevels and therefore does not affect the ground-state 
polarization dynamics. The tensor term of this Hamilto- 
nian is proportional to a rank-two spherical tensor oper- 
ator T2^o ~ 3^z ~ F{F + !)■ A Hamiltonian consisting 
of only rank-one terms acts to rotate the system, while a 
rank-two (or higher) term in the Hamiltonian is generally 
capable of converting between different rank polarization 
moments of the system. 

Consider a single atomic nucleus with a nonzero 
quadrupole moment"^ in a polycrystalline solid. While 
the average electric field "seen" by the nucleus is zero, 
there may be electric field gradients that interact with the 
quadrupole moment according to the single-crystallite 
Hamiltonian (expressed in the Cartesian basis Xi — x, 

X2 ^ y, X3 = z 
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^ In order to possess a quadrupole moment, the nucleus must have 
angular momentum / > 1. In general, a quantum system with 
angular momentum J may have polarization moments of rank k 
ranging from to 2J. 
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FIG. 1: Energy splittings induced in (a) an F = 1 atomic 
system by a uniform electric field, and in (b) an / = 1 nu- 
cleus by the interaction of a quadrupole moment with axially 
symmetric electric field gradients. In both cases, the split- 
ting can result in quantum beats that convert alignment to 
orientation. 



Here Qij is the nuclear quadrupole moment tensor and 
E is the local electric field at the position of the nucleus. 
Upon choosing the Cartesian coordinate system to be the 
principal axis system of the local electric field gradient 
(EFG) tensor, this becomes 



Hq = 
where 



iac^gj [34^ -/(/+!)]+ I (4 +/!)}, (3) 



4/(2/ - 1) dz 



{I,mj=I\Q,,\I,mj^I) (4) 



is the quadrupolar sublevel splitting frequency 
{{I ,mj — I \Qzz\ I ,iTij — I) is the nuclear quadrupole 
moment and dE^/dz is the principal value of the electric 
field gradient tensor) and the quadrupolar asymmetry 
parameter is 



dEx/dx — dEy/dy 
dEjdz ■ 



(5) 



In the following, we assume ry = 0, i.e. that the electric 
field gradients at each nucleus have cylindrical symme- 
try about the z-axis. In this case, the nuclear quadrupole 
Hamiltonian (Eq. I^J is formally analogous to the atomic 
Stark Hamiltonian (Eq. with corresponding similari- 
ties in both the level splitting (Fig. ^ and the dynamics 
of the two systems. 

The interaction (0) lifts the degeneracy between sub- 
levels corresponding to different magnetic quantum num- 
bers |M| of the nucleus. In a sample at thermal equilib- 
rium, the energy splitting gives rise to nuclear polariza- 
tion because, according to the Boltzmann law, there is 
a higher probability of finding a nucleus in a lower en- 
ergy state. ^ Thus the nuclear quadrupolar axis is initially 
along the EFG axis of symmetry. [This is in contrast to 



Typical values of the sublevel frequency splittings are between 
100 kHz and 10 MHz. At room temperature, the relative popu- 
lation difference between sublevels with different |M| is typically 
~10-7. 



FIG. 2: The creation of alignment via (a) optical pumping 
of an _F = 1 atomic system and (b) thermal distribution of 
quadrupole-split 1 = 1 nuclei. In (a), optical pumping with 
linearly polarized light creates atomic alignment along the 
polarization axis of the light. Straight arrows represent light- 
induced excitation and wavy arrows represent spontaneous 
decay. In (b), the alignment axis is determined by the direc- 
tion of the electric field gradients in the crystal. 



the production of atomic alignment by optical pumping, 
in which the alignment axis is determined by the polar- 
ization of the pumping light (Fig. [2J|.] Although each 
nucleus is in an aligned state, in a disordered medium 
such as a powder there is no macroscopic polarization of 
the sample because the distribution of individual crys- 
tallite orientations is random. However, in spite of this, 
NQR signals corresponding to macroscopic magnetiza- 
tion of the whole sample can still be observed in such 
media, as discussed below. 

The initial nuclear alignment of several crystallites 
with different orientations of the local field gradients is 
illustrated in the first column of Fig. |3 using angular 
momentum probability surfaces, as discussed in Ref. 
(see also Ref. Q). For a particular density matrix, the 
distance the surface from the origin in a given direction 
is proportional to the probability of finding the projec- 
tion M = I along that direction. For clarity, we assume 
complete polarization, i.e., that all the nuclei are in the 
lowest energy state. 

The excitation is accomplished by a resonant radiofre- 
quency (rf) magnetic-field pulse. We assume that the 
field B{t) ~ Bi cos(ajt + (/>) (where uj is the rf frequency 
and (j) is the rf phase) has constant amplitude Bi and is 
applied for a time r at an angle (3 to the EFG axis of 
symmetry. We assume that uj is equal to the quadrupo- 
lar frequency ujq, and consider the decomposition of this 
field into components along and perpendicular to the 
EFG axis of symmetry. It can be seen that the longi- 
tudinal component causes rapidly oscillating level shifts 
that have negligible effect on the nuclear polarization, 
whereas the transverse component can be further decom- 
posed into two oppositely polarized circular components, 
each of which drives one of the transitions from M — 
to M' = +1 or —1. Assuming that the Rabi frequency 
LUi — jBi (where 7 is the gyromagnetic ratio) is much less 
than loq, we can neglect the nonresonant component for 
each transition. The resonant components, of amplitude 
Bi sin(/3)/2, cause rotation of the nuclear polarization by 
an angle wiTsin(/3)/2 around the direction of the trans- 
verse component of the magnetic field. This follows from 
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FIG. 3: Probability surfaces 1_7| corresponding to the evo- 
lution of the nuclear polarization in several crystallites with 
different orientations of the (ajcially symmetric) local field gra- 
dients. The angles a and /3 are the Euler angles of the sym- 
metry axes of the local electric field gradients with respect 
to the fixed lab frame. Each row shows the evolution of a 
given crystallite in time. The last row shows an average of 
the polarization over all possible crystallite orientations. The 
plots at times t = Q and r show the nuclear polarizations at 
the beginning and end of the excitation pulse (with magnetic 
field amplitude Bi). After the pulse, a quantum beat cycle 
is shown. As can be seen in the bottom row, macroscopic 
oscillating orientation appears along the direction of Bi . The 
plots are generated from a density matrix calculated using an 
average-Hamiltonian approximation in the quadrupolar inter- 
action frame (see, for example, Ref. [l^'l. The powder average 
is found by integrating analytically over the Euler angles (see 
Appendix). 



consideration of the dynamics in an interaction frame in 
which the quadrupolar interaction is removed and the 
resonant components of the rf field appear to be static 
(see Appendix). Since in a typical NQR experiment the 
pulse length is much longer than the quantum-beat pe- 
riod T = 27r/wQ, quantum beats begin to occur during 
the rf pulse. However, at times when the quantum-beat 
phase is zero, the excitation corresponds to simple rota- 
tion. In the second column of Fig. [SJ we plot only the 
effect of the rotation, and not of the fast quantum-beat 
oscillation, by assuming that the pulse length r is an in- 
teger number of quantum-beat periods. However, none of 
the mechanisms described here depend on this assump- 
tion. We have chosen the parameters of the excitation 
pulse such that the rotation is by 7r/4 for crystallites 
whose EFG axes are orthogonal to the rf polarization 
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FIG. 4: Dynamics of the individual crystallites' alignment 
axes during the excitation pulse. Each dot represents the 
intersection of a crystallite's alignment axis with the unit 
sphere. Initially {uJit = 0) the alignment distribution is 
isotropic (dots are plotted along parallels of latitude and 
meridians of longitude in order to illustrate the dependence of 
rotation on the initial polar angles) . As the phase u}\t accumu- 
lates, the alignment axes rotate (in the average-Hamiltonian 
approximation) by an angle a;itsin(/3)/2 around the direc- 
tion of the component of the rf magnetic field transverse to 
the EFG principal axis, as described in the text. Thus, the 
alignment axes rotate away from and around the rf field axis 
(indicated by the vertical arrow). 



axis.^ 

After the excitation pulse is over, the alignment axes of 
the nuclei have been rotated away from the local electric 
field gradient axes. Thus the nuclei are in coherent super- 
positions of eigenstates of different energies — the requi- 
site condition for quantum beats. These quantum beats 
correspond to a cycle of alignment-to-orientation conver- 
sion, as shown in the last five columns of Fig. |31 In one 
period T — 2'k/ujq of the cycle, alignment is converted 
into orientation, then into alignment at an angle of 7r/2 
with respect to the original alignment, followed by con- 
version to the opposite orientation, and back to the orig- 
inal state. This evolution is the same as the evolution of 
an aligned atomic system in the presence of an electric 
field (see, for example, Ref. 0). 

The powder average over all crystallites is shown in the 
bottom row of Fig. O Initially {t = 0) the sample has no 
average polarization, as indicated by the isotropic proba- 
bility surface. At the end of the excitation pulse of length 
T (where r is chosen to be an integer multiple of T) there 



^ According to common NMR/NQR terminology, the pulse that 
accomplishes such a rotation is called a, -it/2 pulse. The termi- 
nology stems from the two-level spin-1/2 system, where if one 
starts, for example, with a "spin-down" state and applies a pulse 
creating a coherent superposition of "spin-down" and "spin-up" 
states (with equal amplitudes of the two components) this cor- 
responds to rotating the orientation direction by 7r/2. Similarly, 
a TT pulse transfers all atoms from the spin-down state to the 
spin-up state, and rotates the orientation by it. In the present 
case of a spin-one system, if the excitation pulse transfers all of 
the initial M = population into a superposition of the M = ±1 
sublevels, this actually corresponds to a physical rotation of the 
alignment by 7r/2 (not by ttX). Unfortunately, there appears to 
be some confusion in the literature about this point. 
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is a net nuclear alignment within the sample, indicated 
by the elongation of the probability surface along the Bi 
axis. The net alignment arises because the alignment axis 
of each nucleus rotates away from the Bi axis during the 
excitation pulse (Fig. Following the excitation pulse, 
each nucleus undergoes alignment-to-orientation conver- 
sion. Since the orientation produced in each crystallite 
is perpendicular to both the EFG principal axis and the 
axis of the alignment prepared by the excitation pulse, all 
crystallites contribute coherently to the orientation along 
Bi , which leads to a net oscillating orientation® of the en- 
tire sample. This oscillating orientation corresponds to 
a net sample ac magnetization that is the source of the 
detected NQR signal. Details of the calculation used to 
obtain the powder-averaged dynamics, including an an- 
alytic formula for the powder-averaged density matrix, 
are given in the Appendix. 

It must be mentioned that the conversion between po- 
larization moments of quadrupolar nuclei is well under- 
stood in the field of nuclear magnetic resonance in the 
context of multiple-quantum coherences. The NMR sit- 
uation differs from that of NQR in several respects. At 
the high magnetic field strengths common in NMR exper- 
iments, the interaction of the nuclear spin system with 
this field is dominant. Within the rotating frame ap- 
proximation, the quadrupolar interaction for any crys- 
tallite appears to be cylindrically symmetric about the 
external magnetic field axis, and the applied resonant rf 
fields are transverse and appear to be static in the ro- 
tating frame. Additionally, the nuclear spin system is 
initially magnetized (oriented) along the dominant mag- 
netic field direction. When the rf irradiation is weak 
compared to the quadrupolar interaction, the conver- 
sion of orientation to alignment has been recognized as 
the means by which multiple-quantum coherence can be 
created from nuclear spin magnetization during a single 
pulse 0,|23. When strong, short rf pulses are used, dur- 
ing which quadrupolar evolution is negligible, the multi- 
polar (polarization moment) formalism has shown how 
multiple-quantum coherence can be created after a two- 
pulse sequence via orientation-to-alignment conversion 
due to quadrupolar evolution between the pulses plll2^ . 
Furthermore, methods have been introduced to visualize 
the polarization moments of the nuclear spin system in 
terms of graphical representations of the corresponding 
spherical harmonics [23j. 

° Even though each crystaUite is undergoing alignment-to- 
orientation conversion, the sample as a whole has only oscillat- 
ing orientation and no oscillating alignment. Since the system 
is symmetric about the Bi axis, all polarization not along Bi 



Finally, we mention that various techniques for con- 
verting nuclear alignment into orientation have been de- 
veloped for the studies of nuclear moments of short-lived 
nuchdes ^,^2^1. 

In conclusion, we have shown that alignmcnt-to- 
orientation conversion plays a prominent role in the 
phenomenon of nuclear quadrupole resonance, convert- 
ing local nuclear alignment into global orientation, and 
thus causing the appearance of a macroscopic oscillat- 
ing magnetic moment. This has been illustrated using 
the method of angular momentum probability surfaces, 
and the relationship of this mechanism to that of AOC 
in atomic physics has been discussed. In future work, 
it will be interesting to extend the present analysis to 
NQR and NMR in nuclei with / > 1, for which po- 
larization moments higher than alignment are possible 
and the transformations between these moments during 
quantum beats are more complicated than alignmcnt-to- 
orientation conversion |17, 26] . 
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APPENDIX A: THEORY 

The Hamiltonian governing the NQR dynamics in the 
presence of an rf field is derived, for example, in Ref. 
[l8l | . It is convenient to represent this Hamiltonian in the 
{\x), \y), \z)) Cartesian basis^ of the EFG principal axis 
system (oriented at Euler angles a and /3 to the rf field), 
because the nuclear quadrupolar Hamiltonian ||2J) is diag- 
onal in this basis. Transforming to the quadrupolar inter- 
action frame removes the diagonal nuclear quadrupolar 
Hamiltonian and adds additional time-dependent terms 
in the off-diagonal elements. This form of the Hamilto- 
nian, given by 

averages to zero. 

Defined by \x) = {\M = -1) - \M = +1)) \y) = 
-i {\M = -1) + \M = +1)) /V2, \z) = \M = 0>. 



H{t) = TiLOi cos {ut ■ 



2iT? 

cos p e 3 



cos/3eT^"Q* 



Vcosasin/3e"^(^+^)'^«* 



-ism asm 



in/3e-''(i-^)'^«* 



cosasin/3e'(^+S)"«*' 
isinasin/3e*(^-^)"«* 




(Al) 
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facilitates a low-power time-average approximation that removes rapidly oscillating terms. Assuming that the electric 
field gradients are cylindrically symmetric {rj = 0) and that the rf field is resonant with the quadrupolar frequency 
{w = uiq), upon averaging over time-dependent terms the Hamiltonian becomes 

^ / cosasin/3e-"^ ' 

H{t) = —\ isinasin/3e-*"^ I . (A2) 

y cos Q sin /3 e*"^ —i sin a sin f3 e^'^ 

We assume that the initial density matrix (also written in the Cartesian basis) is fully aligned along the EFG axis: 




p(0) = 0. (A3) 



The interaction-frame time dependence is given by 



cos^ Q:sin^(ici;Jt) —icofias\nas\Y?{^uj'^t) — ^ cosasin(aj'j^i) e~"^\ (A4) 
i cos a sin a sin (^w^i) s\v? as\v? (J^^'it) ^ sinasin(ti;it) e"**^ ' 

I cos a sm{u)[t) e^'^ \ sin a sin(a;ii) e'"^ cos^ (^0;^*) 



where = Wi sin /3. In order to find the dynamics of 
the powder average, we transform out of the interaction 
frame and rotate the coordinate system to coincide with 
the lab frame, with Bi along the quantization axis. The 
algebraic form of the resulting density matrix is compli- 
cated. However, since in the powder the only preferred 



direction is along Bi, the powder average will be sym- 
metric about this axis. Thus, if we represent the density 

matrix in the (| — 1), |0), |1)) Zceman basis, the averaged 
density matrix will be diagonal, and we can ignore the 
off-diagonal terms: 



^ — cos oj it — {1 + cos oj[t) cos 2 [3 + Asm u'^t sin fi sin {wt + (j)) 

pLa.h{t) = o I • • • 8 C0S2 (iwit) COS^/? 

3 — COS uj'it — {1 + COS u)[t) COS 2 13 — 4: sin uj[t sin f3 sin {tot + I 



(A5) 



To obtain the result for the powder as a whole, we average over the Euler angles, evaluating the integrals using 
Mathematica software: 



Pavg(i) 



J^"" da phabjt) sin l3 dp 



da sin /? 

'lO + 4:F'{-\Loft'^) -3'Kll_i{uJit)+8F"{-\ujlt^)u)itsin{Ljt + (l)) 0> 

'i + 4F"'{-luft^) 

,0 10 + 4 F' {- {ujlt^) -STTli.iiuJit) -8 F''{- {iolt^) LOit sin {ujt + (, 



Here 



F'{z) = 1F2 



(2);(i^);. 



, F"{z) = iF2 



(2);(y);. 



F"'{z) = 1F2 



(A6) 



(A7) 



pFq[a; b; z] is the generalized hypergeometric function where 



(c)„ = c(c+l)...(c + n-l) (A9) 
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is the Pochhammer symbol, and H„(z) is the Struve func- 
tion, which satisfies the differential equation 

^'y" + ^y' + {^'-n')y = lj^^_^ (aio) 

where 

n\\ = n{n-2){n-A)... (All) 
represents the double factorial. 



This gives the powder average during the rf pulse. The 
terms proportional to sin [uot + 0) represent the quantum 
beats at the quadrupolar splitting frequency. In order to 
find the evolution after the pulse, we set the factors uj\t 
to a constant corresponding to the phase accumulated 
during Rabi oscillation induced by the excitation pulse. 
In the results presented here, we set this constant equal 
to 7r/2. The only remaining time dependence is then in 
the sin {wt + (j)) terms. 
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